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Abstract 

A search for QCD instanton (/) induced events in deep-inelastic 
scattering (DIS) at HERA is presented in the kinematic range of low 
XBj and low . After cutting into three characteristic variables for l- 
induced events yielding a maximum suppression of standard DIS back- 
ground to the 0.1% level while still preserving 10% of the I-induced 
events, 549 data events are found while 363126 (CDM) and 435^22 
(MEPS) standard DIS events are expected. More events than ex- 
pected by the standard DIS Monte Carlo models are found in the 
data. However, the systematic uncertainty between the two different 
models is of the order of the expected signal, so that a discovery of 
instantons can not be claimed. An outlook is given on the prospect to 
search for QCD instanton events using a discriminant based on range 
searching in the kinematical region Q2 > 100 GeV^ where the /-theory 
makes safer predictions and the QCD Monte Carlos are expected to 
better describe the inclusive data. 



1 Introduction 

Instantons [|T], are tunnelling phenomena between topologically different 
vacuum states present in non abelian gauge theories which can not be de- 
scribed by perturbation theory. In the case of QCD, instantons (i) induce 
hard processes violating chirality. Deep-inelastic scattering (DIS) offers a 
unique opportunity to discover processes induced by QCD instantons 
because a sizeable rate is predicted within "instanton-perturbation theory" 

*Max Planck Institut fiir Physik, c/o DESY PHI, Notkestr. 85, 22607 Hamburg, 
koblitz@mail.desy.dc 

^invited talk given at the Ringberg Workshop on HERA Physics on June 19th, 2001 



1 




Figure 1: Sketch of a boson gluon fusion process inducing instanton tran- 
sitions. A virtual photon (with 4-momentum q) emitted by the incoming 
electron fluctuates into a g-g-pair, one of which (with 4-momentum q') fuses 
with a gluon {g) out of the proton while the other {current quark) forms 
a hard jet. The /-subprocess is characterized by the negative 4-momentum 
squared Q'^ = — g'^ of the in-coming quark and x' which is defined in analogy 
to the standard DIS variable XBj. 



1^, ^, ^, and instanton events exhibit a characteristic final state signature 

EE- 

At the HERA machine, which collides 27.5 GeV positrons on 820 GeV 
protons, the predicted cross section is large enough to make an experimental 
observation possible, although the expected signal is still small compared to 
the standard DIS background. Therefore, finding discriminating observables 
based on the expected "fire-ball" like topology of /-induced events produc- 
ing a large number of hadrons and combining these variables to a powerful 
discriminant is the key ingredient for a successful search analysis. 

In contrast to the first pioneering studies which were concentrated on 
relatively few variables characterizing /-induced events W, Tl\, the anal- 



ysis by the HI collaboration |2^, |3^ which is presented here uses dedicated 
variables and optimized cuts. 



2 Instanton-induced processes in DIS 

Instanton processes in DIS are dominantly induced by photon gluon fusion 
as sketched in figure 0, where a photon emitted by the incoming electron 
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splits into a g-g-pair. One of these quarks (g') enters into the /-subprocess. 
The basic reaction in this subprocess is 



where qn (qji) denote right handed quarks (anti-quarks) and g gluons. In 
every /-induced event, one quark anti-quark pair of all nf kinematically ac- 
cessible flavours is produced. Charm and bottom quarks can in principle 
be produced, but their cross section is strongly suppressed. Chirality is vi- 
olated by these events with Ax = 2n/. Anti-Instantons also contribute to 
the cross section in the same manner, but here only left handed quarks are 
in the final state. The quarks are emitted isotropically together with a mean 
of (ug) ^ 0{l/as) ~ 3 gluons and fragment into a densely populated band 
of hadrons with a width of about 1.1 units of 1] [|. They are homogeneously 
distributed in azimuth. This band together with the relatively hard jet origi- 
nating from the current quark form the characteristic final state of /-induced 
events. Since in each /-event a pair of strange quarks is produced, in addition 
an increase in the amount of strange hadrons is expected in the band region. 

The actual number of produced hadrons heavily depends on the squared 
centre of mass energy Wj = {q' + of the /-system. It can be expressed 
in terms of Q'^ and x' by Wf = Q'^(l —x')/x', where Q'^ is the negative four 
momentum vector squared of the quark entering the /-subprocess and X IS cL 
variable describing this process in analogy to the Bjorken scaling variable xbj 
in DIS. The instanton cross section can be safely calculated in "instanton 
perturbation theory" only for large values of Q'^ and x'. Towards lower 
values the cross section calculated perturbatively steeply increases. Matching 
the instanton size and distance distributions in the QCD vacuum, which 
are fundamental properties of instantons, to high quality lattice calculations 
allows to derive lower boundaries Q'^ > 113 GeV^, x' > 0.35 above which the 
perturbative calculations are applicable. 

The most recent calculations of the cross section for /-induced events at 
the HERA collider requires an additional cut on > g'2 ^ Q^y2 
further reduce theoretical uncertainties due to non-planar diagrams, which 
are not present in the calculation. In the kinematic region defined by the 
cuts required by the perturbative ansatz x' > 0.35, > Q'"^ > 113 GeV^ 
and general additional cuts x > 10~^ and 0.1 < y < 0.9 the cross section at 
HERA is crl^l^j^ = 29.2^g;^pb [|, |T^. The error given is only reflecting the 
uncertainties of Aj^. The cut > 113 GeV^ has not yet been included in 
the experimental analysis. 

^The pseudo-rapidity ry is defined as 1] = — logtan(0/2), 6 being the polar angle and 
the proton direction being the positive z-axis. 




flavours 
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Figure 2: Example of a simulated instanton event at = 54 GeV^ 



X 



2.97- 10-3, Q'2 = 136.5 GeV^ and x' 
proton is moving in +2;-direction. 



0.48 in the HI detector. The incoming 



3 QCD Monte Carlo Generators 

The instanton induced processes described above are simulated by the Monte 
Carlo event generator QCDINS which works on top of the HERWIG 



T9| generator and calculates the matrix element of the hard subprocess. The 
defauh implementation of QCDINS is used with x' > 0.35, Q'^ > llSGeVl 
After the hard subprocess has been generated, gluons are emitted according 
to the leading logarithm approximation. The hadronization step describ- 
ing the transition of the outgoing partons to hadrons is done using cluster 
fragmentation . 



A detailed simulation of the hadronic final state of the background DIS 



events is provided by the two Monte Carlo event generators RAPGAP ||2T|] 
and ARIADNE [^]. Both generators incorporate the 0{as) matrix elements. 
The main difference is the treatment of gluon emission where RAPGAP 
implements the parton shower approach ||23| while ARIADNE employs the 
colour dipole model 0]. We will therefore refer to RAPGAP as "MEPS" 
and ARIADNE as "CDM" in the following. 



4 Event Selection 

The data were collected in the year 1997 when the HERA machine collided 
Ee = 27.5 GeV positrons with Ep = 820 GeV protons. The accumulated data 
sample recorded with the HI detector amounted to an integrated luminosity 
of 15.8 pb^^. A detailed description of the HI detector can be found in P3| . 
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A scematic view of a simulated QCDINS event in the HI detector is given 
in figure ^j. For this analysis the important detector components are the 
backward SpaCal calorimeter which is necessary for the measurement of the 
scattered positron and the liquid argon (LAr) calorimeter which together with 
the tracking system measures the particles in the /-band and the current jet. 
Because the positron is only identified using the backward SpaCal calorimeter 
polar angles are limited to 156° < 9^ < 176° which restricts the analysis to 
values of < lOOGeV^. 



5 Finding discriminating variables 

The most characteristic feature of /-induced events is the densely populated 
band of hadrons in the laboratory frame due to the isotropic decay of the 
instanton in its rest system. Therefore the measurement of the energy flow 
of the hadronic final state is crucial to the identification of /-induced events. 
The energy fiow in the detector is reconstructed combining the calorimetric 
energy depositions in the LAr and SpaCal calorimeters with the momenta of 
low momentum tracks {pt < 2 GeV) in the central jet chamber according to 
the procedure described in . Also the second striking feature of /-induced 
events, the current jet, relies on this measurement. 

From a theoretical point of view the most interesting observables char- 
acteristic to /-induced events are x' and Q'^ since they are linked to the 
instanton size and distance which show a characteristic behaviour. There- 
fore reconstructing these kinematic variables of the /-subprocess from the 
hadronic final state is paramount to identifying /-induced events. 

The starting point of the reconstruction of the underlying kinematics 
can either be finding the current jet or reconstructing the expected band 



of hadrons in the detector. Both ways have been intensively studied |27 
and a preference to identifying the band first was given. The first step is 
to reconstruct the mean pseudo-rapidity f] of the band by finding the Et 
weighted mean rj of the calorimeter clusters 

— ' clusters 

and define the band byr^^ G [f/ — l.l,f/ + l.l]. This spread of 1.1 rapidity 
units of the band is expected from basic theoretical considerations. 

The second step is to find the current jet which has roughly an Et of 
5 GeV, as can be seen in figure ^ Using the CONE jet algorithm with a 
cone radius oi R = \/Jjf~+~^ = 0.5 allows to identify the jet in about 70% 
of the cases. The reason for the relatively small cone size is the fact that 
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Figure 3: The shape normahzed distributions of three variables which al- 
low discriminating between standard DIS background and /-induced events. 
Shown are the number of charged particles in the /-band nb, the sphericity 
Sph and the reconstructed virtuality of the quark entering the /-subprocess 
Q'^gc- The lines and the associated arrows denote the cuts applied. 



the correct current jet can also be directed into the /-band and thus a bigger 
cone radius is prone to pick up more energy from the band while the cone size 
is still sufficient to find jets and reconstruct their kinematics reliably when 
they are well separated from the band. If the jet is emitted into the same r] 
and region as the band, then the clusters belonging to the jet are removed 
from the band. After the reconstruction of the jet Q'^ can be calculated and 
assuming the invariant mass of the band to be the mass Wj of the instanton 
allows to derive x' (see figure p. 

The tagged band can now also be exploited to calculate properties that 
should be characteristic to the particles produced by an instanton transition. 
That is a high isotropy of the outgoing particles, a high number of charged 
particles which can be seen in the tracking system of the detector and possibly 
an increased strangeness production. 

For this analysis the following three variables have been chosen to charac- 
terize the /-induced events: The reconstructed virtuality of the quark enter- 
ing the /-subprocess Qreo ^^e number of charged particles in the 77-region of 
the band rib and the sphericity of all particles except those belonging to the 
jet, calculated in their centre of mass frame. These variables are shown in 
figure ^ Because of their potential in discriminating /-induced events from 
background and due to the good description of the data prior to any cuts 
these variables are used in a cut procedure to enhance the /-signal over the 
DIS signal. The distributions are generally described within 5% by the stan- 
dard DIS Monte Carlo programs except for the tails which are only described 
within 10% (for the highest multiphcities up to 20%). 

Three additional observables have been chosen which are used to monitor 
the remaining events after the cuts and to see whether they support the 
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Figure 4: The distributions of three variables which are used to discriminate 
between standard DIS background and /-induced events are shown in plots 
a) — c). They are the number of charged particles in the /-band rib, the 
sphericity Sph and the reconstructed virtuality of the quark entering the /- 
subproccss Q'^^^- Plots d) — f) show the transverse energy in the /-band fi'tjb, 
the isotropy variable Ab and the transverse energy of the supposed current 
jet. 



/-hypothesis after after discriminating against the DIS background. The 
variables are the total transverse energy of the band £'t,b; the isotropy variable 
Aft described below and calculated for the band and the transverse energy of 
the jet /?t,jct- These variables arc described not as well as the first three, but 
still the description differs by about only 10%. 

The isotropy variable of the band A = (i^'m.b — -E'out,b)/-E'm,b is defined 
using the quantities 



but 



E 



mm > \pn 

n Hadr. 



max > \Pn ■ 1 

n Hadr. 



where E^ut will be small for a non-isotropic set of hadrons since then there 

— * 

is a direction i onto which the projections of the hadron momenta Pn is 
small and E^^ will be large for the same set of non-isotropically distributed 
hadrons. Therefore A is approximately unity for a non-isotropic distribution 
and vanishes for isotropic distributions. 
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6 Search Strategy and Results 



In order to reduce the standard DIS background, several different combina- 
tions of cuts into the three observables Q^c, Sph and rib are studied and 
evaluated for their ability to separate /-induced events from the background 
by achieving a high separation power S = eiNs/eois, where e is the fraction of 
events left after the cuts. 125 different scenarios where surveyed which rep- 
resent all possible combinations of five different values at which to cut into 
the three observables. The following cuts are used: Uf, > 5, 6, 7, 8, 9, Sph > 
0.4, 0.5, 0.55, 0.6, 0.65 and for Q'^l^: 95, 100, 105, 110, 115 < Q[l^ < 200 GeVl 
Three different scenarios were singled out according to the following criteria: 

(A) The highest instanton efficiency (ejNs) 

(B) High eiNs at reasonable background reduction (eDis) 

(C) Highest separation power S at eiNs ~ 10% 

The values of the cuts chosen in the three scenarios as well as the number 
of events in the data and the two standard DIS Monte Carlo simulators is 
shown in table |l|. The error on the number of MC events includes systematical 
and statistical uncertainties. The systematic experimental uncertainties are 
of an order of 5-10% and the two available standard DIS Monte Carlo models 
based on leading order QCD matrix elements and implementing different 
treatments of higher order parton emissions agree with each other within 
20%. The number of data events is about 20% larger than predicted by 
the DIS models. Note that this is about the difference of the prediction of 
the two Monte Carlo generators and the background estimation is largely 
unknown in this extreme phase space region. Whether the excess can be 
explained by instanton induced processes or is simply due to an incomplete 
implementation of higher order parton emissions or due to another possible 
deficiency in the description of the DIS background, for example, an incorrect 
treatment of heavy quarks, remains an open question. 

Looking at the three variables after the cuts of scenario C (figure |^) shows 
that the three variables which were used in the cut scenario would support 
the /-hypothesis while especially the distribution of /?t,b and i?t,jet disfavour 
it. In both of these distributions data shows a softer behaviour than the 
DIS Monte Carlo programs while QCDINS predicts a harder spectrum. 
Recent theoretical analysis shows, that this might be due to the lack of 
a cut in > 113 GeV^ which suppresses the contribution of non-planar 
graphs in the /-subprocess and which are not implemented in QCDINS. In 
fact, Schrempp and Ringwald have shown that it is possible to soften the Et 
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Figure 5: The six observables after the cuts to enhance the instanton signal 
of cut scenario C. 



Table 1: The three cuts scenarios A, B and C which are used to reduce the 



background. 



Set 


mario 
Q" [GeV^] 


Cuts 
Sph> 


rih > 


eiNS 


CDM 


NS 

pis 
MEPS 


CDM 


1= Events 
MEPS 


Data 


A 
B 
C 


95-200 
105 - 200 
105 - 200 


0.4 
0.4 
0.5 


5 
7 

8 


32% 
21% 
11% 


35 
56 
86 


34 
52 
71 


2469+i^ 
1005^?^ 
363^1 


1084^1^ 
4351^6 


3000 
1332 
549 



spectra of the jet and the band while not affecting the other four distributions 
by a cut against low x and thus, by exploiting the strong correlation between 
X and in QCDINS and suppressing the contribution of events at low 

It will therefore be definitely interesting to extend the analysis to high 
where the theoretical uncertainties in particular on these two variables 
are smaller. In the following section a possible search strategy is presented. 
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Figure 6: Example of a binary tree storing events 1-7 according to their 
position in the x-y plane. The first event is always the top node of the tree 
which has two daughter nodes for events with smaller and larger x- value. On 
the second level the tree then sorts according to y. 

7 A Prospect of a Search at > 113 GeV^ 

Since the expected cross section of /-induced events in the kinematic region 
{xb > 10~^, 0.1 < y < 0.9, > 113 GeV^) is by approximately two or- 
ders of magnitude lower than the DIS cross section (Jdis ~ 3000 pb, it is 
important to gain the highest possible signal to background ratio. This can 
be achieved by exploiting observables characterising /-induced processes. To 
find these observables a large number of promising event variables have to 
be investigated and the sensitivity to systematic details in the modelling 
of the hadronic final state has to be tested. This requires a sophisticated 
and fast discrimination method to find the appropriate combination of event 
variables. Range Searching allows to define such a powerful discriminant. 

7.1 A Discriminant Based on Range Searching 

Events can be classified as signal or background by estimating their proba- 
bility density p at the point of the event in the event-variable phase space, 
employing a Monte Carlo (MC) generator to sample the densities. In the 
case of neural networks (NN) this is done by fitting the probabihty-density 
with the adjusted weights of the neurons. To circumvent this time consuming 
procedure the density at each point can be directly estimated by counting 
the number of background and signal events in a surrounding box V . Given 
the ratio 

p Pin ^ *nv) 

■ p{DIS) 4DIS{V) 
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the probability of an event to be a signal event is D = + €). Compared 
to NN's this method also has the advantage of not extrapolating into phase 
space regions where there are no sample events available. Thus signals from 
data events outside the region covered by the MC simulation can be avoided. 
Counting the number of events in the vicinity of a certain point is a problem 
known as Range Searching. 

Range Searching algorithms have been developed which allow a search 
time ~ log(n), where n is the number of points that have to be searched |17 . 



An algorithm using binary trees to store the events is employed, suitable 
especially for a large number of events and dimensions (i.e. observables) . 
The MC events are successively filled into the nodes of two binary trees, one 
for the signal events and one for the background events, where the criterion 
by which the decision is taken to descend to the left or right of a node is given 
by the value of one of the event variables (see figure |^). While descending 
the tree this variable cycles through all considered discriminating variables. 
After filling, the position of every event in the tree is given by its coordinates 
in the event variable space. Classification of an event is done by searching 
in the trees for all background and signal events in the box V . This is done 
in the same manner as filling the tree. A more detailed description and the 



properties of this method can be found in [28 



7.2 Results 

Starting with 35 variables based on the hadronic final state the best 12 were 
chosen by calculating the discriminant with all 2-combinations (pairs) of the 
initial variables and taking those variables which provide a high separation 
power S = eiNs/cDis demanding an efficiency for instantons of eiNs = 10%. 
The number of considered variables is further reduced by calculating all 5- 
combinations and selecting those with the highest separation power and a 
small systematic variation of the background. The systematic uncertainty 
was obtained by using four standard DIS-MC simulators which were tuned 
to data on representative hadronic final state quantities, in the range > 
100 GeV^ at HERA [Q. The variables forming the best combination is 
shown in Figure |^. The separation power for eiNs = 10% is S = 126. In 
Figure ^ the shape normalized discriminant D is shown for the /-induced and 
the background events, as well as the distribution for D > 0.9 normalized 
to a luminosity of 100 pb~^ which is comparable to that already collected 
by each of the HERA experiments HI and ZEUS. An event sample can be 
isolated where half of the events are expected to be due to instantons while 
the /-efficiency is still 10%. 

In addition the reconstruction of x' and Q'2 was studied before and after 
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Figure 7: The characteristic event variables providing good instanton sepa- 
ration with small systematic uncertainties. Shown is the reconstructed vir- 
tuahty of the quark entering the /-subprocess Qf^.^, the sphericity of the par- 
ticles not associated to the jet in their rest system, the second Fox- Wolfram 
moment of these particles and the event shape variable -Eout,B which is the 
projection of the particle transverse energy onto the axis, that makes this 
quantity maximal (see [Q). Finally the number of charged kaons in the I- 
Band is shown. 

a cut on D. As can be seen from figure |^, a cut on the discriminant at 
D > 0.95 improves the reconstruction of both variables strongly. The events 
which were selected by the discriminant are very /-like and therefore the 
reconstruction based on ideal /-events is successful. 

8 An Outlook on Proving Chirality Violation 

The one property of /-induced events which is not present in standard deep 
inelastic scattering is chirality violation. Therefore finding events which vio- 
late chirality would be a strong hint for /-induced processes. 

In every /-induced event chirality is violated by Ax = 2n/, meaning that 
all the outgoing quarks are right handed (or all left handed, in the case of 
anti-instantons). In addition, in every /-induced event a pair of s-s quarks is 
produced. The decay of A baryons might provide a unique way of measuring 
the polarization of these strange quarks, provided both quarks form a A 
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Figure 8: To the left, the shape normahzed discriminant D for the instanton 
events using QCDINS and for standard DIS events using four MC simulators 
is shown. The second plot shows a zoom into the rightmost part and is now 
normalized to a luminosity of L = 100 pb~^. At eiNs = 10% 178 /-induced 
events are expected. 



and A-baryons after hadronization respectively. However, the problem here 
is that a large fraction of the A-baryons contain s-quarks stemming from 
fragmentation. So far, no efficient selection mechanism has been found to 
select A-pairs containing only strange quarks produced in the hard instanton 
process. Inferring the polarization of the s-quark from the A-decay has been 



well studied in polarized deep inelastic scattering |31 



An additional problem of such an analysis is the relatively low detector 
efficiency for A-baryons which is of the order of 25% for a rapidity range of 
\ri\ < 1.3 for the Hl-detector, which leads to only a 6% chance of finding 
both A-baryons in an event. 



9 Conclusions 

HI has conducted a search for /-induced events in the kinematic range of 
XBj > 10~^, 0.1 < y < 0.6 and 6^01 > 156° using 1997 data. Using a simple cut 
technique the standard DIS background has been suppressed up to the 0.1% 
level while retaining 10% of the /-induced events. More events were found in 
the data than predicted by the DIS Monte Carlo programs. Nevertheless this 
excess is not conclusive since it is of the order of the difference of the Monte 
Carlo predictions. While four out of the six distributions that were studied 
would support the /-hypothesis after these cuts, the spectra of the band 
and the jet do not favour it. However, recent theoretical analysis has shown 
that this discrepancy can be explained by the lack of a cut > 113 GeV^ 
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Figure 9: The reconstruction of x' and Q'"^ strongly improves when a cut 
on the discriminator is apphed. With a cut D > 0.95 the RMS of the 
reconstruction of x' is 0.20 and 0.17 for Q''^. 



which suppresses non-planar graphs not present in the QCDINS Monte Carlo 
simulator. 

It is therefore of high interest to extend the analysis to this range of to 
reduce the theoretical uncertainties. A first survey of the prospects of such a 
search was done using the Range Searching method to discriminate between 
the /-signal and the DIS background. The results are very promising since 
this method allows to choose combinations of observables with small system- 
atic uncertainties in the final discriminant and a high separation power. 

A possibility of measuring the chirality violation of the /-induced events 
via the decay of A-baryons is in principle possible if an efficient mechanism 
to select A-baryons containing only strange quarks from the hard /-process 
can be found and experimental difficulties can be overcome. In any case such 
an analysis will require a very large luminosity and an isolated data sample 
where the contamination from standard DIS processes is low. 
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